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Abstract 

The iterative algorithm recently proposed by Waxman for solving eigenvalue problems, which 
relies on the method of moments, has been modified to improve its convergence considerably without 
sacrificing its benefits or elegance. The suggested modification is based on methods to calculate 
low-lying eigenpairs of large bounded hermitian operators or matrices. 
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Recently, Waxman [l] has proposed a convergent iterative algorithm for obtaining so- 
lutions of the eigenvalue problem, which does not involve a matrix diagonalization. For 
operators which possess a continuum as well as a set of bound states this is most advan- 
tageous pfl. In the case of the ground state, for example, the eigenenergy, e, is determined 
numerically as a function of the coupling constant of the potential, A, and inverted to yield 
the e corresponding to the required value of A. The convergence rate of the algorithm, there- 
fore, depends on two factors: the number of iterations required to find an eigensolution for 
a particular choice of e and the number of times this must be repeated in order to determine 
theThe iterative algorithm recently proposed by Waxman for solving eigenvalue problems, 
which relies on the method of moments, has been modified to improve its convergence con- 
siderably without sacrificing its value of e which corresponds to the desired value of A. In a 
recent paper we have shown that for many non-singular symmetric potentials which vanish 
asymptotically, a simple analytical relationship between the coupling constant of the po- 
tential and the ground state eigenvalue exists which can be used to reduce the number of 
times A has to be calculated for a given value of e Q|. Here we show, that it is possible also 
to reduce the number of iterations necessary to determine an eigensolution for a particular 
choice of e. Furthermore the existence of a simple analytical relationship between e and A 



can be used to handle problematic situations that are referred to as pseudo convergence in 
other methods |4|] and which can occur in the Waxman algorithm as well. 

In the Waxman algorithm, eigenpairs are determined as functions of the strength of the 
potential in the following manner Here we shall use the abstract Hilbert Space notation 
for the most part and point out their meaning in a one dimensional coordinate space where 
appropriate. Starting from 

(T - XV)\u >= e\u > (1) 



with 



u(x) =< x\u >, lim u(x) = (2) 

x-| — >oo 



where T is the kinetic energy operator (or more generally a suitable hermitian operator for 
the unperturbed system); A > is the strength parameter of the potential (AV > and 
V(x)— > as \x\ — > oo) and the energy eigenvalue, e (with e < 0), is negative and corresponds 
to a bound state. Using Green's method a solution to equation (pQ) is given by 

\u >= XG e V\u > (3) 
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where the Green's operator, G € , can formally be defined as 

G e = (T-e)- 1 



(4) 



and the corresponding Green's function satisfies 

lim G e (x) = 0. (5) 

x-| — >oo 

Normalizing \u > with a suitable reference state \ref > 

< ref\u >= 1 (6) 
allows A to be written as (see equation (J3j) ) 

A =< ref\G e V\u >~ x (7) 
which can then be used to eliminate A from equation (jBJ) 

GeV\u > 

W>= <ref\G e V\u> (8) 

From equations (J2J) and (jSJ), A can be determined as a function of e in the following manner. 
For a particular choice of e equation ([Sj) can be iterated 

G e V\n > 
< ref\G £ V\n > 

until it converges and A can then be determined from equation (j2j). Repeating for differ- 
ent values of e yields a set of different values of the potential strength A. When enough 
points have been determined, a simple interpolation procedure can be used to determine the 
dependence of e on A. 

If we assume, that u is square integrable and therefore can be normalized we can choose 
| it > as a reference vector in (jSJ) yielding 

\u >< u\G e V\u >= G £ V\u > (10) 

from which it is clearly seen, that we are dealing with an eigenvalue problem for the operator 
G £ V with eigenvalue A -1 =< u\G t V\u >. The equation corresponding to j9]) is now given 



by 



G t V\n > . . 

n + 1>= ^77] • 11 

< n\G e V\n > 



From (TlTj) it immediately follows <n|n+l>=l and therefore 

\n + 1 >= \n > +c^\n± > (12) 

where < n\n±_ >= 0, < n±\n± >= 1 and ^ as long as \n > is not an eigenvector. 
Normalizing \n + 1 > yields 

|„ + 1 >'= (1 + {cifp{\n > +c^\n ± >). (13) 

Equation (jT3j) defines the start vector for the next iteration step and it can be easily 
seen, that the relative contribution of \n > is reduced in each iteration step. It is now 
essential to note, that in the n+1 iteration step we have to deal with two vectors: |n > 
and \n± > which define (in general) a two dimensional subspace. Therefore one could ask, 
does equation (TTBl already define the best choice within this subspace to achieve optimal 
convergence. Remember we are looking for a linear combination of the form 

\n + 1 > — ci\n > +c 2 G e V\n >= d n \n > +d^\n ± > (14) 

where the c's and the d's are chosen to satisfy the normalization and the convergence re- 
quirements. Clearly, in the original Waxman scheme C\ = 0. This corresponds to using the 
method of moments or power method [5|, |6| to solve equation (fTO]) . 

Another choice would be to diagonalize G e V projected onto the 2-dimensional subspace 
spanned by |n > and \n± > and take one of the two eigenvectors as the new start vector in the 
next iteration step. Such an approach has been proposed to calculate low-lying eigenvalues 
of an hermitian operator Q and generalized to low lying eigenstates of unbounded hermitian 
operators Q| . 

The iteration scheme is then the following. From G e V\n > calculate a normalized vector 
orthogonal to the vector \n > 

\n± >= {G e V\n > - < n\G e V\n > \n >) * || . . . (15) 

such that < n_i\n± >= 1, < >= and || . . . || is the C 2 norm. The matrix representation 
of G e V in the subspace spanned by |n > and \n± > is 



(16) 
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Figure 1: Typical convergence rates of the coupling constant A for different ground state energies 
e for Waxman algorithm (power) and the modified algorithm (2x2) using a 20x20 matrix model 
Hamiltonian. 

with 

e n =< n\G e V\n >= A^ 1 , v n =< n\G e V\n± >, a n =< n±\G e V\n± > (17) 

Now the 2x2 matrix ( fT6l ) can be diagonalized yielding two eigenvalues and the corre- 
sponding two eigenvectors. One of the eigenvalues always lies above e n the other one below 
e n as long as v n ^ 0. (Here we make use of the fact, that A is real.) If we always chose the 
upper eigenvalue, the e n form a monotonically increasing sequence bounded by the highest 
eigenvalue of G e V. Therefore the sequence is convergent. (Similar arguments hold, if one 
chooses the lower eigenvalue in each step in case the spectrum is bounded from below.) Since 
in each step the modified iteration step is an optimization with respect to the original power 
method step, one might expect an improved convergence rate. 

In order to investigate the convergence properties we have performed a number of calcu- 
lations with matrices where a discrete spectrum for T was chosen and a random potential 
V was added. The signs in equation (pQ) were chosen such that A is positive and the mono- 
tonically increasing iteration scheme (i.e. increasing in terms of e n = A^ 1 and therefore 
decreasing in terms of A n ) described above was chosen. Figure 1 shows a typical conver- 
gence pattern for a 20x20 matrix using the original Waxman method (power) compared to 
the modified scheme proposed here (2x2). The number of iterations needed until a specified 
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Figure 2: Convergence rates of the coupling constant A of different ground state energies e for the 
Waxman algorithm (power) and the modified algorithm (2x2) using a 20x20 matrix which yields 
very poor convergence as model Hamiltonian. 

convergence limit is achieved for the new scheme is roughly one half that required for the 
original Waxman algorithm. Figure 2 illustrates the convergence properties of the modified 
iteration scheme in another, more dramatic, case where the original algorithm has consid- 
erable difficulty to converge. It can be seen from both examples, that the convergence rate 
improves significantly. 

To judge the advantage of the faster convergence one has to take into account that because 
of the last term in equations (ITTjl the new scheme is always one additional iteration step 
(G e V\n± >) ahead and in each iteration step there are more vector operations. While the 
latter are negligible as far as computational ressources are concerned, the additional iteration 
step has to be taken into account. Even then this increases the total number of iterations 
only by one and still results in a considerably accelerated convergence. 

It is interesting to understand why achieving convergence in the second example is so 
slow. In Figure 3 we show the convergence of A n for different energies as a function of 
iteration number n. It can be seen, that convergence in certain cases is delayed by a number 
of steps in a such a way that seems to indicate that e (or A) appears to converge to an 
incorrect value. Such behavior, which we is referred to as pseudoconvergence is common 



in several iterative a 
Lanczos algorithm ]4\ 



gorithms and has been discussed earlier for the Lanczos and modified 
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Figure 3: Convergence curves of the coupling constant A for different ground state energies e for 
Waxman algorithm using a 20x20 matrix as model Hamiltonian which causes psuedoconvergence. 

In general two problems arise with pseudoconvergence: How to detect and how to remedy 
it? The remedy in the present case may be to switch from the power to the 2x2 scheme. 
Since the 2x2 scheme itself suffers from pseudo convergence, further measures as discussed 
in |4| may be appropriate. 

The task of detecting pseudo convergence is more difficult since if you look e.g. at the rate 
of change of the eigenvalue in the successive iteration steps this can only indicate that there 



may be convergence. However, as has 
smooth relations between e and A exist 



3een shown recently for a wide class of potentials 



^|. We therefore show the corresponding relations 
for the pseudo convergence exhibited in Figure 4. It can be seen, that the power method 
exhibits marked deviations from the smooth behavior at e- values 7 and 8. In Figure 2 
we can see, that for these values the power method seemed to have converged rapidly. In 
fact, and we checked this with the exact results, pseudoconvergence occurred which the 
algorithm had failed to detect. Similarly in Figure 4 in the 2x2 scheme one would suspect 
that pseudoconvergence occurred for e = 6 which is actually the case. 

To summarize, the examples suggest that the modified algorithm has considerably im- 
proved convergence rates in general. In the case where pseudoconvergence occurs the savings 
may become dramatic. The use of the Waxman algorithm enables one to detect the occur- 
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Figure 4: Dependence of the numerically determined coupling constant A on the ground state 
energy e for Waxman algorithm (power) and the modified algorithm (2x2) using a 20x20 matrix as 
model Hamiltonian which causes pseudoconvergence. The two curves which are identical when full 
convergence is achieved have been vertically shifted apart from each other for better visibility. 

rence pseudoconvergence reasonably quickly from the marked deviations from the smooth 
dependence of A on e. As in the power method scheme, the modified scheme does not require 
an explicit matrix representation in a large basis. Thus a major advantage of Waxman's 
method, namely the fact that the scheme can be applied directly to operators and wavefunc- 
tions in coordinate space, either using a discretized numerical or a parametrized analytical 
representation jsj], is preserved when migrating to the modified algorithm. 
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